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Abstract 



We study the number of flux vacua for type JIB string theory on an orientifold of the Calabi- 
Yau expressed as a hypersurface in WCP^[1, 1, 2, 2, 6] by evaluating a suitable integral over the 
complex-structure moduli space as per the conjecture of Douglas and Ashok. We show that 
away from the singular conifold locus, one gets a power law, and that the (neighborhood) of 
the conifold locus indeed acts as an attractor in the (complex structure) moduli space. In the 
process, we evaluate the periods near the conifold locus. We also study (non)supersymmetric 
solutions near the conifold locus, and show that supersymmetric solutions near the conifold 
locus do not support fluxes. 

Flux compactification has become an important area of work in string/M-theory enlarging the 
known families of allowed manifolds that would support fluxes, from Calabi-Yau's and tori, to non- 
Calabi-Yau and non-Kahler manifolds [1]. The statistical counting of flux- vacua has been conjectured 
in [2] to be given by a suitable integral over the complex structure moduli space. In this short note, 
we consider type IIB string theory compactified on an orientifold of a compact two-parameter 
Calabi-Yau expressed as a hypersurface in WCP^[1, 1, 2, 2, 6], and evaluate the conjectured integral, 
away from and near the singular conifold locus. Our work, we believe, addresses for the first time, 
the crucial issue whether the moduli space integral so central to the Ashok-Douglas conjecture, 
is solvable in the context of compact Calabi-Yau's whose mirrors involves more than one complex 
strucutre modulus. The same becomes important given the fact that flux compactification in string 
theory can hope to solve the vacuum selection problem by fixing some or all the moduli which in turn 
implies fixing of more than one (around twenty) undetermined parameters of the standard model. 

We show that the integral gives a power law for the number of allowed flux vacua, and that 
the conifold locus acts as an attractor in the complex-structure moduli space, giving the dominant 
contribution to the fiux-vacua counting. Surprisingly, the result for the latter is very similar to 
that of the one-parameter Calabi-Yau expressed as a hypersurface in WCP^[1, 1, 1, 1,4] in [5]. We 
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also investigate the compatibility of supersymmetry with fluxes, and flnd that one can not have 
supersymmtric solutions that also support fluxes near the conifold locus. 

According to the conjecture of Ashok and Douglas, the total number of no-scale flux vacua is 
given by a moduli space (Ai) integral (^^^^)^^ J_\4 det{—R — u), where X4 is the elliptically-fibered 
Calabi-Yau four-fold in F-theory to which the the type IIB orientifold can be uplifted using the 
prescription of [3]; X4 for the case at hand is a Calabi-Yau hypersurface in WCP^[1, 1, 2, 2, 12, 18], 
for which xiX^) = 19, 728[4]. 

Consider the (moduli) space FD{t) x Mcs: where FD(t) implies the fundamental domain of 
the SL{2,Z) group with the complex Teichmuller space parameter r: |i?e(r)| < |, |r| > 1, and 
Aics is the moduli space of complex structure deformations of a Kahler manifold. We begin with 
the evaluation of Jjn^^ det{—R — u), Aics being the complex-structure moduli space with (f),ilj the 
relevant moduli and a; being the Kahler form, for << 1 and around = 2. We show that one gets 
a power law for some positive r. For points near the singular conifold locus: 0= 1 — 864'0^ (one 
can similarly consider = — 1 — 864'^^), the aforementioned moduli-space integral gives a logarithmic 
result (constant) l/(/rir -|- constant) when integrated around |0| << 1, |1 — 864-0^ — 0| << 1. 

The deflning hypersurface for the Calabi-Yau is: 
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12'll)XQXiX2X^Xi — 20X^X2 = 0, 



(1) 



with /i^'^ = 2 and /i^'^ = 128. This is invariant under Z2 x Zg, and using the Greene-Plesser 
construction, under this modding, one gets the mirror manifold with h}'^ — 128 and /i^'^ = 2. It 
is thir mirror that we will be considering, or equivalently, as in [4], the moduli that can appear at 
higher orders, are consistently set to zero. Under the symplectic decomposition of the holomorphic 
three-form VL canonical homology {Aa,B"',a = 1,2,3) and cohomology bases («„,/?"), defining the 
periods as /^^ Q = z", J^a ^ = Fa, such that Q = z^aa — FaP"". Then, the Kahler potential K is given 



ln{—i JcY^ Ail) = ln{—i{T — r)) — ln{—%^YXX), 11 being the six-component 
I3 \ 

In the vicinity of -0 = and on some regular locus, the period vector 11 is given by: 



by: — /n(— i(r — r 
period vector and S 
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In (2), 



which can be analytically continued to ^ \\ dC,^^2j=M=. We will also require 



(<A-C)'^ 



If x = — 2, then using: 

128 VQEllipticK{^) -32^/6 xEllipticK{^) -9 Tlx 2i^i(f , J, 2, |) 



384 TT 



-(-576 ,F,{^,^,l,\) + 48x 2Fi(3^,^,l,i) + 7a; 2Fi(iU,2,i 



^-i(2) 



576 23 



one gets: 



g(-l + (-l)^)7r2(-576 ,Fi(^,^,l,|)+48x ,Fi(^, X, i, i) + 7x ,Fi(f|, ff, 2, 1)) . 
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2^r(|) 

2l((3 + 2i)+4(-l)^ + (-l)A)7rf (-5762^1(3^,^,1, i)+48x2Fi(^,^,l, i)+7x2Fi(i|, if, 2, i)) 

^72""^ r(|f 

-108V'2(-128y6E//ipticif(-) + 32y6a;E//iptic/s:(-) + 9 7ra; 2i^i(^, -, 2, -)))}, 



■ (^-l|)^^(-576 2^1(3^,^,1,!) + 48.x ,Fi(^,^,l,|) + 7x ,Fi(l|, If, 2, |)) , 

2^r(|f 

_ 2I Trf (-576 2Fi(^,^,l,|)+48x 2Fi(3^, ^, 1, i) + 7x 2Fi(ff , f|, 2, |)) 



36 r(|)^ 



-54^^(-128^E//ipticir(^) + 32 V6a;£;ZHpticX(^) + 97ra; 2i^i(^, ^, 2, ^)))}, 



2§7ri(-576 ,Fi(3L,X,i,i) + 48x ,^1(3^, X, 1, i) + 7x ,Fi(l|, 1|, 2, |)) 

^72""^ r(|)' 

54V'^(-128V6£;ZZipiicir(^) +32 V6x£;ZZipticX(^) + 97rx 2i^i(^, ^, 2, ^)))}, 
{Z!;,(2^(l + (-l)'^)^^(-576 2^1(3^, ^,l,i)+ 48a; 3^1(3^, 3^, 1, i) + 7a; 2Fi(i|, if, 2, i)) 



72 r(|f 

^) + 97rx 2^1(^,^,2,1)))}} 

(6) 



2 2 5 T 1 

-108 V'^ (-128 £;iiipiicis:(-) + 32 x EllipticK {-) + 97rx 2Fi{-,-, 2, -)))}} 



Using the numerical values of various quantities, the expression of the Kahler potential excluding 
the axion-dilaton modulus, is given by: 



K — —ln{—i{T — t) — In 



(1.65278 X 10*^ + 2.32831 x 10"'" i) + (762417 - 546469 i) V - 



(163034 + 2.91038 x 10"^^ i)x- (236957 - 169840 i)i)'^x + (762417 + 546469 i) + 
-(75206.6 + 53904.9 i) a; - (2.32831 x 10"^° i) |^|^ a; + -(163034 - 2.91038 x IQ-^U) x 
-(75206.6 - 53904.9 i) + (16082.1 - 1.81899 x 10"^^ i) |a;p + (23373.9 - 16753.4 i)^^ |^|2 



4 



-(236957 + 169840 i)ip'^x + (2.32831 x 10"^° i) x + (23373.9 + 16753.4 i) V' 



1 2 72 



= —ln{—i{T — t) + In A+ {aiip^ + a\ip^) + {h^x + feix) + {ciip^x + ciip'^x) + {c2Xip^ + C2X'0^ 



+(c3|x| V +C3|a:;| ■0 + C4|a;| 



(7) 



where we drop Im{A) ~ 10 /m(c4) ~ 10 and a CslV'l^x + C5|'0|^x term with Rec^ — 0, Imc^ 
10~^^^. From (7), one gets the metric: 
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2V'(t-#) 

V 2V^(f-^) + ^ / 



(8) 



One then can calculate the curvature 2-form: R^^^^dz^ A dz'' using that for a Kahler manifold, 



dV 



^I'ki ^ ^ = l^ya'^^i.djgin + digjn))- So, using the metric of (8), one gets: 
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^Of course in^SII] is real - the fact that one gets infinitesimal, but non-zero imaginary parts as well, seems to 
be an artefact of using Mathematica. 
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.jf. — djd]^det{g,fj), is a hermitian matrix, i.e., Rjfj = Rjf. or equivalently 



is a hermitian form. However, 7?™ 



R"},dz^ A dz'' is not a hermitian form.^ 



One can show that R 
iddlndet{gfj 

For the Kahler metric corresponding to the Kahler potential of (7), one getsii?^^. = djdiln 
i.e., 



R 



'jk 




(11) 



is hermitian, as expected. 

One has to evaluate det{R + cu) which means one has to evaluate the determinant of: 
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(12) 



One thus gets: 



/ det{R + u;)=[ {R^^^ + grf)dT A df 

JfdxMcs Jfd 



A 



L 



+{r: 



Mcs 



FD 

V-V-V ^ 9^^)i^''xxx + 9xx) - {R\xrP + 9-il>x){R''x^a: + 9xi,) " (-^tci/iV ^ 9xi,){R'',j,xx + 9-ipx) 



xxip + 9xx){R 



- -I- R'^ F?^ - 4- R^ _ R^ 

^ VV'a: a;a;i/' iixx^ xipip ^ xipx i'xip 



^We thank D.Joyce for discussfon on this point. 
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— - 



dip A dip A dx A dx 



(13) 



One thus sees that moduh space integral factorizes into an integral over the fundamental domain 
of S'L(2, Z), and an integral over the complex structure moduli space. Given that the former will 
give only a multiplicative contribution, in the following, we will concentrate only on the complex 
structure moduli space integral. 

Assuming that one performs the angular integrals first and then the radial integrals, one gets: 



/ det{R + uj)^ [ 

J J\x 



{\ilj\')\x\d\x\d9MdW0^ 



4„2 ^x—J\p-L 



(14) 



It is interesting that the power is exactly equal to 2/i^'^(CY3) + 2. 

For points near the singular conifold locus: 0=1 — 86Aijj^, the period vector 11 is given Y: 
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where WiS are to be determined as follows. Using [6, 4], in the neighborhood of the conifold locus: 
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1 6 
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(16) 



where 



n=l 



r[6(n + i)] 



_ ~ (r[n + i])3r[3(n+|)] 

r[6(n+|)] ^^'^^ 



-0). 



(17) 



7 



. Now using the Stirling asymptotic series for the gamma function: r{z) 
C^^^^, and the following expressions for Ui,{(j)) for large |i/|: 



\ z \ e 



1 + Tk + 



Uu{(t>) 



Ifc"!— +00 



;i + (/.r+^-i(0-ir+5 



Uu{-(t>) 



>oo 



(18) 



one gets the following asymptotic expression for 

d'^Wi 9 q(1 - (t>) 



dijj^ 7r2-02 (1 - 864V'6 - 0) 
Coj = — > stn — ]stn{ — 



, where 



i(2j — l)7vr 



_ \ ^ • 2 

C2j+1 = V 2^ 
r-=l 



vrr \ f Tir\ i(2j+i)7rr 



■(-r- 



One can show that q = (1, 1, —1, —2, 2, 1)[6, 4]. ^ One thus gets: 



Ci /27ri (1 - 864^'^ - 0) 
2^V4^ (1 - 0)2 



Ml- 864^^6 -0) + /,(0,V^), 



(19) 



(20) 



where fi{(f>,ip) are analytic functions of and ■0- Now, the monodromy properties, as discussed in 
[7, 4] imply that: 

Wi = ^(wo - wi)ln{l - 864^^^ - 0) + /,(0, V). (21) 

This imphes Wq — Wi — fo — fi- To determine we use the fact that for|0|<<l, |'0|<<1, one 
can expand fi as: fi = J2n,m (^n,m(/^""4^^ [7]- In the neighborhood of the conifold locus, the analytic 

part of the periods /i(0, -0) will be given by: fo{e^^^'^'^, e^t/j), j — 0, 5, where /o(0, "0) is given by 
{x = l- 864V^6 _ 



do,: 
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^If one includes the in the Stirhng expression for V{z)^ one wih get a /n(l — 864^/;^ — 0) term in addition to the 
1 in ^ 

(1-864j/'<'-0) ■ 
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In general, 



1^ r(f)(-i2V')"^-^(0) 



/o(0,^) 6£r(n)(r[i-f]mi-i|]' 

where 



Now, from (20), one sees that fo — fi = ax + Pxcf) + x = 1 — S6At/j^ — 0, and 
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The numerical values of the relevant infinite series are: 



^2V "°-"^ /„" ^-11.6-0.5i, 
n=l (864)6 

^2g aon(-)tn ^ 

n=i (864)^ 

12V '\"°'V^ ^-13.3-1.4i, 
=1 (864)? 



n= 

OO 



^l?. (SM)! + 

^^v^ ea ao (-). ^ 
„=i (864)6 

12 y — = -12.1 + 24.4i, 

n=l (864)6 

12 V = -34.2 - 25.9Z, 

n=i (864)^ 

OO i2i2i / \ n 
■5-^ e 2 fln « — 6 n 

"g (^f ^-82.5 + 7.7i. 

12 y , °'" y — = -7.1 - 82.5i, 

(864)6 



9 



Similarly, 
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From the asymptotic analysis done earlier, writing fi = af^^ + aj^V + (li^^x + (9(0^, x^, one 
knows that /q - /i = {a'^^ - af^)x + 0{(j)'^,x'^,x(l)). From (28), one sees that this is approximately 
satisfied if one drops the imaginary part of /o (to an accuracy of about 95%) and also the imaginary 
part of /i (to an accuracy of 89.5%) and assumes that 0— >0,a;— >0:^— >0. One then gets, 
ft) — fi = -^ix + C(0^,x^,0a;) (up to 94.5% accuracy). Given that the periods can be regarded 
as the solutions to the homogeneous Picard-Fuchs equation, this implies that a constant times the 
period vectors will also be valid period vectors. Consequently, one can do a further rescaling of the 
periods: Wi —>■ to get a match with (20). 

The Kahler potential is thus given by: 

K = -ln(-i(r - f) - ln(-m^J:U) = 

= (^l/o - fi\M^\' + Re(^[3/o + 3/i + /2 + /a - /4 - h][-fo + /2 + A] 



—ln(—i{T — r)) — In 

+ [/0 + /2][/0 + /l-/2 + /3-/4 + /5]) 




f ) + a\x\Hn\x\' + A + {Bx + Bx) + {C<p + C<p) + D\xY + 

G02), (29) 

where S = ( V ^^"^ have considered terms in the Kahler potential up to quadratic in 

X, (and their complex conjugates and the products of the same). Using gfj — didjK and from our 
experience from points away from the conifold locus, we realize that we can exclude the axion-dilaton 
modulus as it will have only an overall multiplicative effect in the moduli space integral. Thus, one 
gets: 




(30) 



Writing 



cs 01 (\CY . ,2 (\CY E\(\BY {a^D)\ \CB\^ 



one gets: 



^^'^ =-iIa^-IJ^^I^I +V^^ A)\A- A ) A- 

= Ai/n|a;|' + A2, (31) 



rdx A dx R'^^^A,dx A dx 
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1 / (^^-^^^dxAdx 

(Ai/n|x|2 + A2)2 1 _c%Ai 



(32) 



which even though is not hermitian, the corresponding R^j = didjlndetg = didjln{ ^^^^^.^t^^^ ) , hence 
given by: 

f ^ (T-f)2 
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(33) 



is hermitian. Dropping the Kahler form as compared to the curvature two-form, as x — >^ 0, 

Constant r|(Constant) 



f det(R + uj)^ / \x\d\x\\(l)\d\<P\-——— j—^--^ (34) 

One thus sees, that the contribution from the conifold locus to the complex structure moduli space 
integral dominates over the contribution from points away from the singular conifold locus, as ex- 
pected. The conifold locus, thus, acts as an attractor in the complex structure moduli space. 
Lets now consider the (non)supersymmetric solutions, implying that we look for solutions to: 

W = DrW = D^W = D^W = (35) 

W being the superpotential, given by: 

[ {F-TH)An = {2'Kfa\f - rh) ■ H, (36) 

where one uses that the NS—NS flux H and the RR flux F are given by: F — (27r)^Q;'(J\j/3a-|-J\j+3Q;a) 
and H — {27rya' {TiaPa + Ti-a+scta), cta, P"" forming an integral cohomology basis, a — 1,2,3. For 
nonsupersymmetric solutions, W ^ 0. Now, W = together with 

DrW = drW + drKW = --^—{J^^U - tH^U) = 0, (37) 

(r -r) 

implies 

JT^n = Ti^U = 0. (38) 

The analysis below is applicable for nonsupersymmetric solutions, however for supersymmetric so- 
lutions, after doing the analysis below from (41)- (43), one notes that if in and H^, one sets 
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J-i = Ha = 0, then one can write U = 7i U = These conditions translate to: 

-^^1 + Y(3/0 + 3/i + /2 + /3-/4-/5)+^3(/0 + /2) + y(-/0 + /2 + /4) 
+ Y(/0 + /i-/2 + /3-/4 + /5)=0 

-^^1 + ^(3/o + 3/i + /2 + /a - /4 - h) + Hsifo + /s) + ^(-/o + /2 + A) 

+ ^(/o + /l-/2 + /3-/4 + /5)=0, 



which imphes 



(39) 





+ af 


+ 4°' 


-4°^-4°^) + ^3( 


;4n4°^) + f( 


-4°^ + ar 










-4°^ + 4°^) = o 








^(3a« + 3aW 


+ a?^ 




-4')-4')) + ^3i 


;4n4^^) + f( 










+ 4^^ 


-4'^ + 4'^) = o 









+ ^(3a?^ + 3af ) + 4'^ + 4'^ - 4'^ - 4'^) + ^3(4'^ + 4'^) + $(-4'^ + 4'^ + 4'^) 



2^ , 2 \ — u ' 1 ' "Z ' "O "O / ' - o\"U ' "Z / ' 2 

+ Y(4'^ + 4^^ - a?' + - af + 4^^) = 0. (40) 

Similar equations can be written for 2,3,5,6- It is understood that equations (41)-(43) are to solved 
for x,(f) in terms of JF, 7i and r. One has to look for integer solutions (for jFj and Tii) of (40). For 
instance, to solve for jFj, one sees that one ends up with three equations in five variables. Keeping 
any two fixed, one in general gets complex solutions in terms of the these fixed TiS. One thus has 
to ensure that the imaginary parts vanish and the real parts turn out to be integral. One can show 
that it is not possible to achieve the same. This implcs that near the conifold locus of the (orientifold 
of) the compact Calabi-Yau expressed as a dcgrcc-12 hypcrsurfacc in WCP"'[1, 1, 2, 2, 6], one can 
not obtain supersymmetric solutions that support fluxes. Note, however, for non supersymmetric 
solutions, for H'^U ^ 0, JT^H ^ 0, D^W = implies that r = 
We now consider D-^VF = 0. 

D^W = d^W + d^KW = 
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27r 



,(2) 



.(2) 



a 



(2) 



(2) , „(2) 



,(2) 



+ 1) + 4" 

I 



On 



+(47r(xZn|x|^ + x) + S + + F0)(:r^ - rTi^) 



' 2-K 



\ 



U^f0 + Sfi + f2 + fs-f4-h) 
/0 + /2 

l^/nx + /o 

i(-/0 + /2 + /4) 
V i(/0 + /l-/2 + /3-/4 + /5) y 



+ (3 a^o^ + 3 ai^^ + a^^^ + a^' 



(2) , ^(2) 



(2) 



(-a?^ + + a':>) (^5 - ^5 r) + (4^^ + a'. 



>(2) 



(2) 



.(2) 



a^'^) (^2 - r) + 2 (4'^ + ) (J^s - ^3 r) + 



(6 4 



(1) 



,(1) 



,(1) 



An 
.(1) 



/nx , , —i {T\ — Hi t) 



TT 



(2) 

03 X — 



4 X — as X + 6 ag + 02 + 4 ~ '^4 ~ '^5 0) (-^2 — ^^2 t) + 



4" 



Ig , 
.(1) 



X + a^^^ X + 4'^ + a^'^ 0) (J^3 - ^) + 



,(1) 



(1) 



a^fP X + 0,2^'' ,x + a)^' X — ag'"'' + 02'"'' + 04'"'' 0) — Tis r) 



,(2) 



,(1) 



(1) 



,(1) 



+(24- 
,(1) 



(1) 



4'^ 



,(1) 



.(2) 



I (2) 



(2) , (2) (2) , (2) , r, (1) 

a2 X + ay X — al X + X + 2 Qq ' ' 



a"^'' + ai'^ - 4'^ + al'' 0) (J^g - Tig r) 



,(1) 



(1) 



+2 (:r4 - t) (4'^ + o?^ 2; + 4'^ + ^=)) (S + 2//,x + F0 + x(L' + 47r + 47r /nlxH) ) = 

4 TT y 

Near x — 0, = 0, and for nonsupersymmetric solutions (for which ^4 0, Ha ^ 0) one gets 
4 7r2 / ^ (^1 - 7i 1 r) (3 ^ + 3 af ^ + a^^ + - - af ) (^2 - ^2 r ) 



,(1) 



xlnx 



(41) 



^4 — 7^4 T 



TT 



a, 



(2) j_ „(2) 




,(2) 



4') +4^0 (-^5- ^5 r) 



,(2) 



+ 4^0 (-^3-^3 r)- 
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.(2) 



+ ai 



(2) 



af+ (-5^6 -Her) 



.(2) 



B -i {Ti - Hi z) 



+ (6 + + - - af) (^2 - r) + 



,(0) 



,(0) 



TT 



(0) 



,(0)^ 



2 (4°) + 4°)) (^3 - ^3 r) + 2 (^4 - 7i4 r) + (-4°) + af + af) (^5 - ^5 r) + 



(42) 



For nonsupersymmetric solutions, for which ^4 = 7^4 = 0, one gets a hnear equation in and x from 
(41); there will be no Inx or Inx terms. 
Lets now consider DjW = 0: 



/ 

i(34^) + 34'^ + 4'^ 



.(1) 



4^4^^ 

a? 



i(-4'^ + 4^ + 4 



(1) 



2 1^0 



+(C + E(j) + Fx + 2G(j)){T^ - tH^) 



1(3/0 + + /2 + /3 - /4 - /5) 



■^Inx + /o 

K-/0 + /2 + /4) 



V \{h + h-f2 + h-h + h) j 

^ ((3 4') + 3 a? + a? + 4'^ - 4'^ - 4'^) {^2 - 7^2 r) + 2 (4^^ + 4'^) (^3 - r) + 



2 4'^ (^4 - ^4 r) + (-4^^ + a'^'> + al'>) {T^ -n^T) + {a\;> + a[ 

{C + 2Gcl> + Fx + e^) ^ -'^^^-^^^) + (64°) + 4°) + af - af^ - af^ + 



(1) ^ ^(1)^ 



,(1) 



.(1) 



- 4') + 4^0 (-^6 - ^6 



TT 



3 Cq^^ X + 3 X + X + 03^^ X — af^ x — a'^^ x + 6 a^'^ (j) + 



a^^ + al'' - - a^'^ 0) , (^^3 _ r) + 2 (a^" + a^" + 4'^ x + a^'^ x + 



,(1) 



,(1) 



,(1) 



(0 



,(0 



,(2) 



4'^ + 4'' 0) (^3-W3r) + 
4'^ 0) (^5-H5r) + (2 4°)- 



,(1) 



.(0) 



,(2) 



(2) 



,(0) 



+ Oo^-* x + a^^'' a; — a'2' x + 03^'' x — a^^' x + 05^'' x + 2 g'q' — 



,(2) 



4°^ 



4°^ + 



(2) 



,(2) 



X + a2 X + a\ X 



(j) + a^2^ + 
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a^2^ 4) + ai^ (t> - a^P + a!^^ 0) (J^e - Tie r) + 2 (J^4 - ^4 r) (ai°^ + x + + 



xlnx 




47r2 



= 



Substituting (42) in (43), one gets a single complex constraint on and (j). 

The fact that we are able to show that the conifold locus acts like an attractor in the complex 
structure moduli space is what can also be argued qualitatively given that the moduli space integral 
involves the curvature tensor and hence the contribution to the same will be dominated from locii 
where the curvature tensor diverges. As part of future work, one could try to solve equations (41) - 
(43) numerically and perform a Monte-Carlo simulation like [5] to verify, e.g., that the conifold locus 
indeed involves clustering of solutions, and hence acts Uke an attractor in the moduli space.. 

After completion of this work, we became aware of [8], which has overlap with our work. 
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